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INTRODUCTION 


Several  years  ago  the  Naval  Weapons  Center  (NWC)  embarked  on  an 
effort  to  establish  expertise  In  the  areas  of  radoae  analysis,  design, 
and  testing.  An  interdepartmental  program  was  devised  to  develop 
structural,  thermal,  and  electromagnetic  capabilities  at  NWC  and  to  up¬ 
grade  existing  test  facilities.  As  a  result  of  this  program,  the  RF 
and  Microwave  Technology  Branch,  Code  3814,  has  become  knowledgeable  in 
the  area  of  electromagnetic  radome  analysis  and  has  created  a  number  of 
computer  codes  which  can  predict  many  useful  properties  of  radomes. 
Starting  with  the  simplest  of  models  and  progressing  toward  more  com¬ 
plex  variations,  we  have  built  up  a  large  library  of  subroutines  which 
can  be  combined  in  various  ways  to  give  a  very  flexible  radoae  analysis 
capability. 

This  report  contains  a  description  of  the  general  theory  for  ana¬ 
lyzing  radomes.  The  various  types  of  computer  codes  are  described  in 
the  section  on  Theory  and  Discussion.  Examples  of  the  Inputs  and  out¬ 
puts  of  each  code  are  presented  in  the  section  which  carries  that 
title.  The  Summary  discusses  all  progress  to  date,  and  the  expectation 
of  work  to  be  done  is  detailed  under  Future  Work. 

At  the  time  of  this  writing,  our  electromagnetic  analysis  effort 
will  continue  for  at  least  one  more  year.  Thus,  all  work  is  ongoing 
and  subject  to  change. 


THEORY  AND  DISCUSSION 


The  electromagnetic  radoms  codes  can  be  divided  into  four  general 
groups.  Each  code  In  a  group  Is  completely  self-sufficient,  yet  is 
related  to  the  codes  In  all  other  groups.  The  first  group  predicts 
boreslght  error  as  a  function  of  antenna  gimbal  angle.  The  second 
group  computes  reflection  and  transmission  coefficients  for  a  multi¬ 
layer  system  of  flat  (no  curvature)  surfaces.  The  third  group  computes 
various  parameters  characterising  two-dimensional  Von  Karmen  shaped 
radomes,  and  the  fourth  group  consists  of  codes  based  on  the  propaga¬ 
tion  of  spherical  waves  through  spherical  radomes* 
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HUKKSIGHT  ERROR  CODES 

The  boroslght  error  codes  calculate  boreslght  error  ns  a  func¬ 
tion  of  antenna  glmbal  (look)  angle  using  the  technique  of  plane  wave 
ray  tracing.  There  are  three  codes  In  tills  group:  (1)  two-dimensional 
horeslght  error  transmit  formulation;  (2)  two-dimensional  borcslght 
error  receive  formulation  using  monopulsc  processing;  and  (3)  three- 
dimensional  borcslght  error  receive  formulation,  also  using  monopulsc 
processing. 

The  transmit  case  is  based  on  the  following  theory.  An  antenna 
is  placed  inside  the  radome  and  divided  Into  equidistant  point  sources 
(Figure  l),  each  of  which  emits  a  ray  traveling  perpendicular  to  the 
antenna.  The  points  at  which  each  ray  intersects  the  radome  are  com¬ 
puted,  the  normals  at  these  points  are  found,  and  the  incidence  angles 
at  which  each  ray  strikes  the  radome  are  calculated.  The  flat  panel 
transmission  and  reflection  coefficients  are  used  for  multiple  iso¬ 
tropic  layers,  and  the  raya  are  refracted  as  they  propagate  through 
varying  dielectric  constants,  which  differ  from  that  of  free  space. 
The  rays  are  then  extended  to  a  constant  phase  plane  outside  the  radome 
where  a  best-fit,  least-squares  technique  is  used  to  form  a  straight 
line  through  the  points  In  this  phase  plane.  The  difference  between 
the  tilt  angle  of  this  line  and  the  glmbal  angle  of  the  antenna  is  the 
bo resight  error. 

The  receive  formulations  (both  two-  and  three-dimensional)  are 
very  similar  excqpt  that  an  Infinite  plane  wave  Is  assumed  incident 
from  outside  the  radome  (see  Figure  2)  and  Is  divided  Into  equidis¬ 
tant  points.  Each  point  emits  a  ray  traveling  perpendicular  to  the 
Incident  plane  wave  which  Intersects  the  radome.  As  in  the  transmit 
formulation,  the  intersection  points,  normals  to  the  radome  at  those 
points,  incidence  angles,  and  transmission  and  reflection  coefficients 
arc*  computed.  The  rays  are  then  extended  to  the  antenna  which  is 
divided  in  half  for  two  dimensional  and  in  quadrants  for  three  dimen¬ 
sional.  The  received  signal  amplitude  is  computed  over  each  half¬ 
quadrant  of  the  antenna,  and  the  monopulsc  equations  are  used  to  com¬ 
pute  the  boreslght  error.  The  two-dimensional  formulation  provides 
only  principal  plane  boresight  error  prediction,  while  the  three- 
dimensional  receive  formulation  in  Figure  3  provides  information  for 
all  planes  (see  Appendix  A). 


PLAT  PANEL  CODES 

Because  the  flat  panel  transmission  and  reflection  properties  for 
multiple-layer  radomes  are  such  an  important  part  of  the  previous  bore- 
sight  error  codes,  we  have  crested  two  programs  which  calculate  these 
cumulative  complex  transmission  and  reflection  coefficients  as  func¬ 
tions  of  frequency  with  incidence  angle  as  a  parameter.  One  code 
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assumes  thnt  all  layers  are  Isotropic  (see  Figure  4),  while  the  second 
m<  Hu  met*  thnt  all  layers  are  uniaxial  ly  anisotropic  (Reference  1),  l.e. , 
the  dielectric  constant  along  the  axis  normal  to  the  surface  Is  differ¬ 
ent  from  that  In  the  plane  of  the  surface  (see  Figure  5).  Both  codes 
t  Irst  calculate  the  single-surface  reflection  coefficients  for  both 
parallel  and  perpendicular  polarizations.  Then  the  flat-panel  coeffi¬ 
cients  for  each  layer  are  computed  (Reference  2).  The  layers  are  then 
cascaded  using  st„  al  flow  graph  theory  to  provide  the  overall  complex 
transmission  and  reflection  coefficients  for  multiple  layers  (see 
Appendixes  B  and  C). 


VON  KARMAN  GEOMETRY  CODES 

The  majority  of  our  electromagnetic  analysis  to  date  has  concen¬ 
trated  on  Von  Karman  shaped  radomes,  although  the  codes  are  set  up  In 
such  a  way  that  extensions  to  any  desired  shape  are  easily  Imple¬ 
mented.  As  a  result,  there  are  two  codes  which  calculate  various 
values  characterizing  a  two-dimensional  Von  Karman  radome.  One  com¬ 
putes  the  functional  value,  first  derivative  or  slope,  second  deriva¬ 
tive,  and  radius  of  curvature  as  functions  of  position  along  the  radome 
(see  Appendix  D).  The  other  computes  the  radome  radius  and  the 
clearance  of  an  antenna  that'  is  gimballed  inside  the  radome  at  any 
point  along  the  axis.  This  code  can  handle  not  only  a  constant  wall 
thickness  between  the  inner  and  outer  shapes  but  also  a  linearly 
tapered  wall. 


SPHERICAL  WAVE/SPHERICAL  RADOME  CODES 

Finally,  In  our  efforts  to  progress  toward  more  complex,  physi¬ 
cally  satisfying  models,  we*  have  developed  two  programs  based  on 
propagation  of  spherical  waves  from  a  dipole  antenna.  One  uses  an 
approximate  spherical  wave  ray  tracing  method  (Reference  3)  and  calcu¬ 
lates  the  electric  field  magnitude  ns  a  function  of  polar  angle  for  a 
dipole  offset  along  the  axis  of  a  spherical  radome.  The  aecond  code 
calculates  exactly  the  same  quantities  but  uses  the  boundary-value 
method  to  compute  the  exact  solution  (Reference  4).  These  two  solu¬ 
tions  are  then  compared  to  show  those  conditions  which  limit  the 
accuracy  of  the  spherical  wave  ray  tracing  approach  (Reference  5). 
While  these  two  codes  are  applicable  only  to  spherical  radoMs*  their 
comparison  has  produced  insights  into  the  limitations  of  spherical  wave 
ray  tracing  which  can  be  applied  to  other  more  typical  redone  shapes* 
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TNPUTS  AND  OUTPUTS 


HO RE SIGHT  ERROR  COOKS 

The  two-dimensional  boresight  error  codes  require  Inputs  specify¬ 
ing  the  radotne  shape,  material  parameters  characterizing  the  various 
layers,  antenna  size  and  frequency,  and  polarization  information. 
Specifically,  the  user  must  input  frequency,  length,  and  diameter  of 
the  Von  Karman  shape,  antenna  position  along  the  radome  axis,  antenna 
diameter,  dielectric  constant,  loss  tangent  and  thickness  for  each 
layer,  total  number  of  layers,  number  of  points  into  which  the  antenna 
length  is  divided,  and  parallel  or  perpendicular  polarization.  The 
three-dimensional  code  requires  all  of  the  above  Inputs  with  roll 
angle,  antenna  polarization,  and  incident  plane  wave  polarization  as 
additional  parameters. 

Examples  of  the  output  of  these  three  codes  are  given  in  Plgures  6 
through  8,  while  comparison  with  recent  experimental  data  taken  in  our 
radome  test  chamber  is  shown  in  Figure  9.  Reasonable  correlation 
between  theory  and  experiment  at  X-band  frequv'.  ;s  has  been  achieved 
using  the  plane  wave  ray  tracing  approach,  althuugh  certain  discrepan¬ 
cies  have  been  found  at  lower  frequencies.  These  discrepancies  are 
being  investigated  currently.  Figures  10  and  11  are  outputs  from  the 
boresight  error  transmit  formulation  code.  Average  and  central  ray 
transmission  magnitudes  versus  gimbal  angle  are  plotted  from  a  modified 
version. 


FLAT  PANEL  CODES 

The  flat  panel  codes  which  plot  reflection  and  transmission  versus 
frequency  require  inputs  consisting  ot  minimum  and  maximum  frequency 
values,  a  frequency  increment,  angle  of  incidence  (in  degrees),  number 
of  layers,  and  dielectric  constants,  loss  tangents,  and  thicknesses  of 
each  layer. 

Examples  of  the  output  of  the  isotropic  code  are  given  in  Figures 
12  through  14.  *  These  codes  provide  an  exact  (no  approximations) 
analysis  of  any  multilayer  isotropic  or  unlaxlally  anisotropic  flat 
panel  system. 


■k 

Examples  of  outputs  from  the  unlaxlally  anisotropic  flat  panel 
code  are  not  included  since  they  would  look  virtually  identical  except 
for  dielectric  constant. 
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SPHERICAL  WAVE/SPHERICAL  RADOME  CODES 

Inputs  to  the  spherical  radome  programs  are  inner  and  outer  radii, 
dipole  offset  distance,  wavelength,  index  of  refraction  (real  only), 
and  observation  radius.  These  codes  plot  the  electric  field  magnitude 
as  found  by  both  the  ray  tracing  method  and  the  exact  method  as  a  func¬ 
tion  of  the  polar  angle  6.  Comparison  of  the  two  techniques  is  shown 
in  Figure  15,  assuming  increasing  source  offset  distance.  Comparison 
of  the  exact  theory  with  experiment  is  shown  in  Figure  16  (Refer¬ 
ence  4). 


SUMMARY 


An  overview  describing  computer  codes  for  electromagnetic  analysis 
and  design  of  radomes  has  been  presented.  These  codes  provide  a  broad 
base  for  predicting  the  electromagnetic  properties  of  radomes  without 
requiring  costly  and  time-consuming  testing  procedures.  As  such  they 
serve  as  an  efficient  design  tool  for  establishing  an  Initial  baseline 
radoroe  configuration.  We  have  created  nine  different  codes  which  com¬ 
pute  boresight  error,  reflection  and  transmission  properties  based  on 
flat  panel  models.  Von  Karman  geometry  characteristics,  and  spherical 
wave  propagation  effects.  Plane  wave  ray  tracing  has  been  the  standard 
approximation  technique  upon  which  the  boresight  errci  codes  are 
based.  The  remaining  codes,  except  for  the  spherical  wave  ray  tracing 
code,  are  all  exact  solutions  which  can  be  used  in  many  different  types 
of  problems.  We  have  tried  to  develop  a  versatile  and  flexible  radome 
analysis  capability  which  can  be  used  in  practical  applications. 


FUTURE  WORK 


Many  relatively  small  modifications  are  planned  in  the  next  phase 
of  the  radome  analysis  and  design  program.  We  will  modify  the  bore¬ 
sight  error  codes  to  include  weighting  functions  which  reflect  differ¬ 
ent  antenna  aperture  distributions.  Both  amplitude  and  phase  weighting 
will  be  considered.  Different  radome  shapes,  curvature,  multiple 
reflections  (ghost  lobes),  complex  permeabilities,  tapered  radomes,  and 
nonconstant  wall  thicknesses  will  be  included  eventually  in  the  bore¬ 
sight  error  codes.  Elliptical  polarization  will  be  built  into  the  flat 
panel  codes.  The  effects  of  tuning  or  taping  a  radome  to  reduce  its 
boresight  error  will  be  analyzed. 

Several  discrepancies  among  the  boresight  error  codes  have  been 
discovered  at  frequencies  below  X-band,  i.e.,  2  to  8  gigahertz.  We  are 
currently  examining  our  monopulse  processing  methods  and  deriving  more 
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rigorous  equations  for  this  aspect  of  the  two-  and  three-dimensional 
receive  formulations  which  should  correct  this  problem. 

Finally,  the  spherical  wave /spherical  radome  codes  will  be 
modified  to  include  multilayers,  and  we  will  use  the  spherical  wave  ray 
tracing  code  as  a  starting  point  for  our  own  more  rigorous  Huyghen's 
principle  spherical  wave  ray  tracing  method  (see  Figure  17). 
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FIGURE  5.  Multilayer  Uniaxially  Anisotropic  Flat  Panel  System. 
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FIGURE  13.  Differencial  Phase  Delay  Versus  Frequency 
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INNER  RAD  =  11.81  SPHERICAL  RADOME 

OUTER  RAD  =  12.70 

INDEX  IS  1.500 

DIPOLE  OFFSET  =  0.200 

OBSERVATION  RAD  -  13.00 

WAVELENGTH  =  1.181 


ELECTRIC  FIELD  MAGNITUDE 


FIGURE  15.  Effect  of  Increasing  Dipole  Offset  Distance. 
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WAVELENGTH  =  1.18 


INNER  RAD  -  6.00 
OUTER  RAD  =  6.89 
INDEX  IS  1.26 

DIPOLE  OFFSET  *  0.94  -  COMPUTED  EXACT  SOLUTION 

OBSERVATION  RAD  =  7.00  - EXPERIMENTAL  EXACT  SOLUTION 


FIGURE  16.  Comparison  of  Experimental  and  Computed  Exact 
Solutions. 
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Appendix  A 

COORDINATE  SYSTEM  FOR  THREE-DIMENSIONAL  RADOME 

In  radome  coordinates,  we  have 


k  *  propagation  vector 
k  -  projection  of  k  in  the  xy  plane 

*y 

k  »  projection  of  k  in  the  xz  plane 

k  »  projection  of  k  in  the  yz  plane 

6R  “  roll  angle  from  y-axis  to  k^ 

By  *  yaw  angle  from  kxz  to  z-axis 
Bp  “  pitch  angle  from  kyZ  to  z-axis. 


Using  direction  cosines, 

A  A  A  A 

k  ■  cosa  x  +  cos  8  y  +  cosy  z 
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where  y,  the  look  angle,  and  6r,  the  roll  angle,  are  given,  the 
relationships  between  direction  cosines  and  pitch,  roll,  and  yaw  angles 
are  given  by 

cosa  *  k  »  k  sing  -  singsing 
x  xz  y  y 

(A-2a) 

or 

cosa  ■  k  sing  *  sinysing 
xy  K  K 

(A-2b) 

cosg  -  k  -  k  cosg_  -  sinycosgD 
y  xy  R  R 

(A-2c) 

or 

cosg  *  k^singp  ■  sinasingp 

(A-2d) 

cosy  ”  k  -  k  cosg  “  singcosg 
'  z  xz  y  y 

(A-2e) 

or 

cosy  »  k  cosg-  «  sinacosg- 
yz  r  r 

(A-2f) 

Thus,  solving  for  the  pitch,  roll,  and  yaw  angles 
cosines,  we  have 

in  term  of  direction 

6^  -  tan-1 (cosa/ cosy) 

(A-3a) 

gR  -  tan-1 (cos a/ cosg) 

(A-3b) 

g  -  tan- 1 ( cos g/ cosy) 

(A-3c) 

28 


Alternatively,  we  could  write 

cosa  -  [cot20R  +  cot20^  +  l}-1^2 

cosg  -  [tan2BR  +  cot20p  +  l]-1^2 


(A-4a) 

(A-4b) 


cosy  “  [tan20y  +  tan20p  +  l]-1^2 


(A-4c) 


The  matrix  which  transforms  coordinates  from  the  antenna  (double- 
primed)  system  to  the  radome  (unprimed)  system  is  given  by* 


-sin0  cos0 
y 

sin0 

-COS0  CO80 

y 


sin0  sin 0 

y 

COS0 

cos0y8in0 


x  " 


(A-5) 


We  have  assumed  a  two-axis  glabal  system  where  Equation  A-5  is 
derived  by  rotating  about  the  original  Y-axis  followed  by  a  rotation 
about  the  new  X-axis. 
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Appendix  B 

MULTILAYER  ISOTROPIC  FLAT  PANEL  EQUATIONS 


Using  Figure  4  which  shows  6  as  the  angle  of  incidence  for  a 
plane  wave  striking  a  plane  sheet  of  Material  composed  of  different 
dielectrics  of  different  thicknesses. 


cos6  -  /Ci(l-j  tanfi^)  -  sin20 
cos0  +  /e^l-j  tan«  )  -  sin20 


perpendicular 

polarization 


I 


c.(l-j  tan6  )co80  -  /e  (1-j  tanfi  )  -  sin29 
— - 1 - - - - -  (B-2) 

E^l-j  tanfi^ )cos0  +  /e  (1-J  tanfi^)  -  sin20 
parallel  polarization  ,  i  »  1,2,3. ..N 


The  rj  (for  both  polarizations)  are  called  the  single-surface 
reflection  coefficients  or,  alternatively,  the  Fresnel  formulas. 
These  form  the  reflection  and  transmission  coefficients  for  a  single 
layer,  l.e.. 


Ti  “ 


(1  -  r2)i 


-J*idi 


1  -  r 


-J28ldi 


(B-3) 
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I 


(B-5) 


and  A0  is  the  free-space  wavelength. 

For  a  aingle  layer,  these  coefficients  can  be  cascaded  using 
signal  flow  graph  theory  to  fora  overall  reflection  and  transaission 
coefficients  for  aany  layers. 

For  two  layers,  we  have 


as  the  signal  flow  graph  representation  which  can  be  rewritten  as 
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(B-6) 


(B-7) 


(B-8) 


(B-9) 


(B-10) 


CB-11) 
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Appendix  C 

MULTILAYER  UNI AX I ALLY  ANISOTROPIC  FLAT  PANEL  EQUATIONS 


Using  Figure  5  for  a  plane  wave  Incident  on  a  multilayer  sheet 
of  unlaxially  anisotropic  material,  we  have 


ri  “ 


^iz  7  811,20  "  co80/elxelz 


parallel 


(C-l) 


/e.  -  sln28  +  cosd/e,  e. 

lx  lx  lz 


cos 8  -  /eix  -  sln28 


cose  +  /cix  -  sin20 


perpendicular 

1  m  1,2, ...N 


(C-2) 


These  single-surface  reflection  coefficients  In  turn  form  the 
reflection  and  transmission  coefficients  for  a  single  layer.  They  are 


(1  -  r|)i 


"Jkldl 


(C-3) 


R. 


-J2k  d 

^(1  -  e  1  *) 


1  -  r*e 


■*J  2kf  di 


(C-4) 
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where 


k  /e  (c.  -  ainz6) 
o  ixz  1Z 


parallel 


(C-5) 


k  Jt, -  sln26  perpendicular  (C-6) 

o  ix 

i  -  1,2, ...N 

with  kQ  «  2n/X0  and  eixz  “  €lx/Eiz*  &l  an<®  can  then 

be  used  to  form  the  overall  reflection  and  transmission  coefficients 
for  many  layers  as  in  Appendix  B. 
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Appendix  D 
VON  KARMAN  GEOMETRY 

In  two  dimensions,  the  Von  Karman  equation  is 


y  -  —  J cos”" 1  ( 1  -  -  j  sinf 2co8-1(l  -  l|l/2  (D-l) 

2/ir  *  * 


where  D  is  the  diameter,  i  is  the  length,  and  h  is  the  antenna  gimbal 
position.  The  first  derivative  or  slope  is 


£ - 5 -  {1  -  cos f 2cos“l ( 1  -  ^±i>]} 


X  { cos" 1  ( 1  -  )  -  -J  sin[2cos-1(l  -  2-(x-£h)->  ]  j"1 /z  (D-2) 


The  second  derivative  is  given  by 

t2v[(-tu(Jl  -  2(x+h))  +  4s  -  u2v2]  (D-3) 

dx*  4/ir 


where 


t  -  {(x+h)U  -  (x+h)J}“l/2 
v  -  {cos-^l  -  \  •in{2co»“l(l  - 


* 
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u  -  1  -  cos(2co8”,(1  -  JiSlill)] 

l 

8  -  8in[2cOS-1(l  -  ] 

l 


and  the  curvature  at  any  point  on  the  Von  Karinan  curve  (Equation  D-l) 
is 


(D-4) 


The  radius  of  curvature  at  any  point  Is  given  by 


R 


K 


{ 
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